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Abstract 

A  long-standing  technological  challenge  to  the  widespread  application  of  thermoelectric  generators  is  obtaining  high-performance 
thermoelectric  materials  from  abundant  elements.  Intensive  study  on  PbTe  alloys  has  resulted  in  a  high  figure  of  merit  for  the  single¬ 
phase  ternary  PbTe-PbSe  system  through  band  structure  engineering,  and  the  low  thermal  conductivity  achieved  due  to  nanostructuring 
leads  to  high  thermoelectric  performance  for  ternary  PbTe-PbS  compounds.  Recently,  the  single-phase  p-type  quaternary  PbTe-PbSe- 
PbS  alloys  have  been  shown  to  provide  thermoelectric  performance  superior  to  the  binary  and  ternary  lead  chalcogenides.  This  occurs  via 
tuning  of  the  band  structure  and  from  an  extraordinary  low  thermal  conductivity  resulting  from  high-contrast  atomic  mass  solute  atoms. 
Here,  we  present  the  thermoelectric  efficiency  of  nanostructured  p-type  quaternary  PbTe-PbSe-PbS  composites  and  compare  the  results 
with  corresponding  single-phase  quaternary  lead  chalcogenide  alloys.  We  demonstrate  that  the  very  low  lattice  thermal  conductivity 
achieved  is  attributed  to  phonon  scattering  at  high-contrast  atomic  mass  solute  atoms  rather  than  from  the  contribution  of  secondary 
phases.  This  results  in  a  thermoelectric  efficiency  of  ~  1.4  over  a  wide  temperature  range  (650-850  K)  in  a  p-type  quaternary  (PbTe)0.65 
(PbSe)0.i(PbS)o.25  composite  that  is  lower  than  that  of  single-phase  (PbTe)o.85(PbSe)0.i(PbS)0.o5  alloy  without  secondary  phases. 

©  2014  Acta  Materialia  Inc.  Published  by  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  world’s  demand  for  energy,  as  well  as  concerns 
about  global  warming,  have  inspired  intensive  research  to 
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develop  clean  and  sustainable  energy  sources.  Thermoelec¬ 
tric  devices  that  convert  heat  to  electricity  are  considered 
promising  candidates  for  waste  heat  recovery  [1],  and  the 
search  continues  to  provide  high-performance  thermoelec¬ 
tric  materials  from  abundant  earth  elements.  A  large  num¬ 
ber  of  recent  studies  [2-8]  have  focused  on  the  mid-range 
temperature  (600-900  K)  thermoelectric  materials,  specifi¬ 
cally  PbTe,  which  has  been  proven  to  possess  a  high  energy 
conversion  efficiency  of  ~1.4  as  defined  by  the  thermoelec¬ 
tric  figure  of  merit,  zT  =  S2To/{ke  +  kl),  where  S  is  the 
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Seebeck  coefficient,  a  is  the  electrical  conductivity,  T  is  the 
absolute  temperature,  and  kl  and  ke  are  the  lattice  and 
electronic  thermal  conductivity,  respectively.  However,  Te 
is  rare  in  the  Earth’s  crust,  in  demand  for  other  applica¬ 
tions,  and  therefore  entails  high  material  costs  for  scale- 
up.  The  focus  of  recent  research  is  to  offer  inexpensive 
alternatives  by  replacing  Te  with  Se  or  S  to  form  binary 
PbSe  [9-11]  or  PbS  [12,13]  or  ternary  systems  of  PbTe- 
PbSe  [4,14,15],  PbTe-PbS  [13,16-18]  and  PbSe-PbS  [19]. 
Nevertheless,  it  is  still  a  big  challenge  to  find  a  Te-reduced 
Pb-chalcogenide  compound  that  can  provide  a  high  zT 
over  a  wide  temperature  range. 

In  a  recent  report  [20],  we  demonstrated  enhanced  ther¬ 
moelectric  performance  of  single-phase  quaternary 
(PbTe)fi_x_>,)(PbSe)J,(PbS)Y  alloys  compared  to  binary 
PbTe,  PbSe,  PbS,  and  ternary  (PbTc)(i_x)(PbSe)x, 
(PbTe)f i_Y)(PbS)Y  systems.  The  electronic  band  structure 
is  tuned  by  soluble  PbS,  the  solubility  limit  of  which  is 
increased  in  the  PbTe  matrix  in  the  presence  of  PbSe. 
The  lattice  thermal  conductivity  is  also  reduced  by  phonon 
scattering  on  high-contrast  atomic  mass  solute  atoms.  A 
more  extensive  report  [21]  on  a  similar  system  (PbTe)(i_2x) 
(PbSe)x(PbS)x>  confirmed  our  findings  for  single-phase 
p-type  quaternary  lead  chalcogenides  and  showed  a  figure 
of  merit  ~2  for  a  single-phase  alloy  with  optimized  compo¬ 
sition  and  carrier  concentrations. 

Recent  studies  [13,16,18]  have  also  shown  that  nano¬ 
structuring  in  the  ternary  PbTe-PbS  system  by  introduc¬ 
ing  secondary  phases  in  the  matrix  is  a  viable  approach 
to  enhance  the  thermoelectric  efficiency.  In  the  present 
work,  we  have  fabricated  bulk  quaternary  Pb-chalcogen¬ 
ide  nanocomposites  through  PbS  alloying  beyond  its  sol¬ 
ubility  limit  in  the  (PbTe)0.9(PbSe)0.i  matrix.  The 
composition  range  of  the  quaternary  (PbTe)(i_x_rf 
(PbScyPbS)x  system  is  large  and  the  choice  of  optimum 
dopant  concentration  will  add  to  this  complication.  Here, 
the  initial  investigation  on  the  thermoelectric  perfor¬ 
mance  of  quaternary  Pb-chalcogenide  composites  was 
limited  to  1  at.%  dopant  concentration  (Na)  which  corre¬ 
sponds  to  Pb0.98Na0.02Te(1_x_>,)Se>,Sx  in  the  equivalent 
formula,  and  the  PbSe  concentration  was  restricted  to 
10  at.%  0  =  0.1)  in  the  (PbTe^i^.^PbSe^PbS)*  sys¬ 
tem  for  all  samples,  to  eliminate  variations  in  the  elec¬ 
tronic  band  structure  of  PbTe  which  can  be  made  by 
Se  solute  atoms  similar  to  the  ternary  (PbTe)(i_x)(PbSe)x 
system  [4,22],  The  thermoelectric  performance  of  com¬ 
posite  samples  is  compared  to  the  parent  solid  solution 
ternary  compound  of  (PbTe)0.9(PbSe)0.i.  The  electronic 
band  structure  alteration  is  only  associated  with  the  frac¬ 
tion  of  PbS  in  the  (PbTe)(0  9_x)(PbSc)0.i(PbS)x  system. 

We  provide  insight  into  the  fundamental  issues  regard¬ 
ing  the  electronic  transport  properties  of  nanostructured 
composite  materials  with  varying  secondary  phase  solubil¬ 
ity  as  a  function  of  temperature.  The  increased  concentra¬ 
tion  of  sulfur  in  the  matrix  at  high  temperatures  and 
variation  in  secondary  phase  fraction  affects  the  electronic 
transport  properties  of  composites. 


2.  Experimental 

2.1.  Sample  fabrication 

2.1.1.  Synthesis 

Polycrystalline  samples  of  PbS,  PbSe  and  PbTe  were 
prepared  by  mixing  high-purity  Pb  (99.999%),  Te 
(99.999%),  Se  (99.999%)  and  dried  S  (99.9%)  in  vacuum- 
sealed  quartz  ampoules  at  a  residual  pressure  of  ~10-4  Torr. 
These  were  reacted  at  high  temperature  to  produce 
high-purity  PbSe  and  PbS  starting  materials.  The  final 
polycrystalline  (PbTe)(0.9_X)(PbSe)0.i(PbS)x  (x  =  0,  0.05, 
0.1,  0.15,  0.2,  0.25)  samples  were  synthesized  by  mixing 
stoichiometric  quantities  of  high-purity  PbS,  PbSe,  Pb 
and  Te,  with  1  mol.%  Na  added  as  the  dopant.  A  total 
mass  of  10  g  was  sealed  in  carbon-coated  quartz  tubes 
under  vacuum,  and  then  heated  to  1373  K  with  a  heating 
rate  of  100  Kir1.  After  being  held  at  1373  K  for  10  h, 
the  samples  were  quenched  in  cold  water,  followed  by 
annealing  at  773  K  for  48  h. 

2.1.2.  Sintering 

The  resulting  ingots  from  the  synthesis  procedure  were 
hand-ground  to  powder  with  a  mortar  and  pestle.  They 
were  sintered  at  773  K  for  1  h  in  a  12  mm  diameter  graphite 
mould,  at  an  axial  pressure  of  40  MPa,  which  was  achieved 
by  induction  hot  pressing  under  an  argon  atmosphere  [23], 

2.2.  Transport  properties  measurements 

2.2.1.  Resistivity  and  Hall  measurements 

Samples  were  loaded  onto  a  heated  BN  substrate  and 
four  probes  were  attached  to  the  edge  of  the  sample.  The 
sample  was  placed  in  vacuum  with  a  magnetic  field  (up 
to  ±2T)  perpendicular  to  its  surface.  The  resistivity  (p) 
and  Hall  coefficient  (RH)  (along  the  hot-pressing  direction) 
were  measured  using  the  van  de  Pauw  method  [24], 

2.2.2.  Seebeck  coefficient  measurements 

The  Seebeck  coefficients  were  obtained  along  the  sam¬ 
ple’s  hot-pressing  direction.  The  samples  were  placed  in 
contact  with  a  heater  on  each  surface  in  a  vacuum  chamber 
[25],  Two  Nb-Chromel  thermocouples  were  also  pressed 
against  the  two  surfaces  of  the  sample  by  spring  force. 
The  heaters  were  programmed  to  provide  a  temperature 
difference  oscillation  of  ±7  °C,  whilst  maintaining  a  set 
average  temperature.  The  thermoelectric  voltage  and  tem¬ 
perature  on  each  surface  were  recorded,  with  the  slope  giv¬ 
ing  the  Seebeck  coefficient  at  the  average  temperature. 

2.2.3.  Thermal  conductivity  measurements 

The  thermal  conductivity  (zc)  was  calculated  from 
k  =  pDTCp.  The  laser  flash  method  (Netzsch  LFA457) 
was  used  to  measure  the  thermal  diffusivity  (Z)r),  the  den¬ 
sity  ( p )  was  calculated  using  the  measured  weight  and 
dimensions.  The  specific  heat  capacity  (C' )  was  estimated 
by  {Cp(kB  per  atom)  =  3.07  +  4.7  x  10“*  x  (T/K  -  300)} 
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[26,27],  which  is  believed  to  be  accurate  for  lead  chalcoge- 
nides  [4,10],  The  combined  uncertainty  for  all  measure¬ 
ments  involved  in  zT  determination  is  ~20%. 

2.3.  Materials  characterization 

2.3.1.  X-ray  diffraction 

The  crystallographic  structure  and  composition  were 
characterized  by  X-ray  diffraction  (XRD)  using  a  PANa- 
lytical  X’Pert  Pro  X-ray  diffractometer  with  Cu  Ka  radia¬ 
tion  (2  =  1.544  A,  40  kV,  30  mA).  In  order  to  measure 
the  phase  ratio  and  to  calculate  the  lattice  parameters, 
the  XRD  patterns  were  refined  using  Rietveld  analysis. 

2.3.2.  Electron  microscopy  analyses 

Scanning  electron  microscopy  (SEM)  was  performed  on 
finely  polished  samples  using  a  JEOL  JSM-7500FA  field 
emission  gun  scanning  electron  microscope.  Samples  were 
also  characterized  by  transmission  electron  microscopy 
(TEM)  using  a  JEOL  2010  microscope.  TEM  samples  were 
prepared  by  cutting  them  into  3  mm  diameter  discs  using  a 
Leica  TXP  polisher,  then  grinding  and  polishing  them  to 
<100  pm  in  thickness.  The  samples  were  then  dimpled, 
and  finally  Ar-ion  milled  on  a  stage  cooled  with  liquid 
nitrogen.  During  ion  milling  low  voltages  and  currents 
were  used  to  reduce  damage  to  the  samples. 

3.  Results  and  discussion 

Samples  with  the  composition  (PbTe)(o  9_X)(PbSe)o  i 
(PbS),,  (x  =  0,  0.05,  0.1,  0.15,  0.2,  0.25)  doped  with 
1  at.%  Na  on  Pb  sites,  which  correspond  to  an  equivalent 
formula  of  Pbo.98Nao.o2Te(o.9_X)Se0.  i  Sx,  were  synthesized 
in  ingots  by  melting  and  annealing,  followed  by  hot  press¬ 
ing  of  the  hand-ground  powders.  The  selected  composi¬ 
tions  are  located  at  the  PbTe-rich  side  of  the  PbTc  PbS 
system  [28],  where  the  phase  separation  of  NaCl-structured 
PbTe-rich  and  PbS-rich  phases  occurs  through  the  nucle- 
ation  and  growth  process.  The  relevant  phase  diagram  is 
discussed  in  detail  in  the  Supporting  Information. 

The  purity  and  crystal  structure  of  the  samples  were 
determined  by  indexing  the  powder  XRD  patterns  in 
Fig.  1 .  The  XRD  patterns  of  samples  with  *  <  0. 1  appear 
to  be  single  phase,  whereas  samples  with  x  >  0. 1  reveal 
two  distinct  phases.  The  inset  in  Fig.  1  shows  that  high- 
angle  diffraction  peaks  are  shifted  with  increasing  PbS  con¬ 
centration  to  10  at.%  due  to  the  alloying  effect.  However, 
additions  above  10  at.%  produced  no  further  shift.  Rietveld 
refinement  was  employed  to  determine  the  lattice  parame¬ 
ters  of  the  matrix  accurately  by  extrapolating  from  high- 
angle  diffraction  peaks.  The  lattice  parameters  of  the 
matrix  and  precipitates,  and  the  proportions  of  secondary 
phase,  are  summarized  in  Table  1.  The  lattice  parameter 
of  the  secondary  phase  is  larger  than  for  pure  PbS 
(a  =  5.93  A),  which  indicates  that  the  PbS-rich  secondary 
phase  is  alloyed  with  PbSe  (a  =  6.13  A)  and  PbTe 
( a  =  6.46  A),  both  of  which  phases  possess  larger  lattice 


10  (degree) 


Fig.  1.  Room  temperature  X-ray  diffraction  patterns  for  (PbTe)(o.9  x) 
(PbSe)0.i(PbS)x  (x  =  0,  0.05,  0.1,  0.15,  0.2  and  0.25)  alloys.  The  inset 
shows  shifting  of  diffraction  peaks  to  higher  angles  for  x  up  to  0.1  due  to 
the  alloying  effect. 


Table  1 

Lattice  parameters  and  percentage  of  precipitates  in  the  matrix. 


Sample 

Lattice  parameter  of 
matrix  PbTe(Se)  (A) 

Lattice  parameter  of 
PbS(Se),  precipitates  (A) 

wt.%  of 
precipitates 

x  =  0 

6.423 

x  =  0.05 

6.406 

- 

- 

x  =  0.1 

6.379 

- 

- 

x  =  0.15 

6.381 

5.945 

4.7  ±  0.2 

x  =  0.20 

6.384 

5.968 

11.1  ±0.3 

x  =  0.25 

6.389 

5.975 

16.8  ±0.4 

parameters.  The  lattice  parameter  of  the  matrix  varies  with 
PbS  addition  up  to  x  =  0.1,  but  remains  almost  constant 
for  x  >  0.1,  as  shown  in  Table  1.  This  indicates  that  for 
samples  with  x>0.1,  the  chemical  compositions  of  the 
matrix  and  the  second  phase  are  fixed  at  any  given  temper¬ 
ature,  although  the  weight  fraction  of  the  second  phase 
does  increase  with  x  above  this  threshold. 

The  SEM  image  of  the  sintered  sample  SI 5  presented  in 
Fig.  2a  illustrates  the  morphology  and  wide  size  distribu¬ 
tion  of  micrometer-sized  precipitates  in  the  matrix.  The 
precipitation  in  multigrain  solids  occurs  preferentially  at 
grain  boundaries  due  to  lower  chemical  driving  force  for 
nucleation  and  higher  diffusion  coefficient  of  elements  at 
grain  boundaries  than  the  grains.  The  segregation  of  solute 
atoms  to  grain  boundaries  may  also  occur  during  solidifica¬ 
tion  [29],  Therefore,  the  precipitates  at  grain  boundaries 
are  generally  larger  than  those  formed  within  the  crystal. 
This  sample  was  exposed  to  several  heating/cooling  cycles 
during  measurement  of  transport  properties,  and  dissolu¬ 
tion  and  coarsening  of  precipitates  may  have  occurred. 
The  corresponding  TEM  micrograph  of  the  same  sample 
in  Fig.  2b  shows  that  an  extensive  population  of  submicron 
precipitates  throughout  the  matrix  is  also  present.  The 
morphology  of  the  precipitates  is  similar  to  the  PbTe- 
PbS  system  [30].  Based  on  the  information  obtained  from 
the  pseudo  binary  phase  diagram  of  PbTe-PbS  [31],  the 
secondary  phases  are  believed  to  be  PbSe-alloyed  PbS. 
Precipitates  within  the  PbTe-rich  matrix  (Fig.  2b)  were 
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Fig.  2.  (a)  SEM  image  of  sintered  (PbTe)o.7s(PbSe)o.i(PbS)o.i5  (S15)  sample,  showing  the  morphology  and  the  distribution  of  precipitates  in  the  matrix;  the 
sample  has  been  through  several  cycles  of  heating  up  to  850  K  and  cooling  during  transport  property  measurements,  (b)  Bright-held  TEM  micrograph  of 
sulfur-rich  precipitates  in  the  tellurium-rich  matrix  when  the  matrix  was  tilted  to  the  [1 1 0]  zone  axis  in  the  (PbTe)o.75(PbSe)o.i(PbS)o.i5  (S15)  ample  shown 
in  (a).  Strong  moire  patterns  are  observed  due  to  the  double  diffraction,  (c)  The  experimental  selected-area  diffraction  pattern  (SADP)  obtained  from  the 
precipitate  ringed  in  (b),  The  diffraction  spots  are  marked  for  precipitates  and  matrix  by  P  and  M,  respectively,  (d)  EDS  spectra  obtained  from  the  matrix 
and  precipitate.  Precipitate  intensities  have  been  normalized  to  match  matrix  intensities  at  the  Pb  L „  line  at  10.55  keV.  The  precipitate  contains  no  Te  and 
is  clearly  a  distinct  phase.  The  Pb  M ^  and  S  Kx  lines  overlap  at  2.3  keV  and  this  peak  is  marginally  more  intense  for  the  precipitate  (Pb  +  S)  compared  with 
the  matrix  (Pb  only). 


characterized  using  electron  diffraction.  Fig.  2c  demon¬ 
strates  the  diffraction  pattern  from  the  [1 1 0]  zone  axis  of 
the  area  circled  in  Fig.  2b,  selected  by  an  aperture  which 
included  matrix  and  a  single  precipitate.  The  precipitate 
has  the  same  crystal  structure  and  orientation  as  the  matrix 
but  its  lattice  parameter  is  smaller.  The  large  number  of 
additional  reflections  in  the  [1 1 0]  pattern  (Fig.  2c)  arises 
from  double  diffraction.  Typically,  single  scattering  of  the 
primary  electron  beam  occurs.  However,  scattered  beams 
can  undergo  further  scattering,  giving  rises  to  multiplicity 
of  the  original  diffraction  pattern,  offset  by  the  original  dif¬ 
fraction  vector.  Fig.  2b  illustrates  the  bright-field  image 
from  which  the  [1 1 0]  diffraction  pattern  (Fig.  2c)  was 
obtained.  The  precipitates  show  very  strong  moire  fringes, 
caused  by  interference  between  beams  scattered  from  the 
matrix  and  precipitate.  All  precipitates  within  a  single  grain 
exhibit  the  same  contrast/fringe  orientation,  indicating 
they  all  have  the  same  orientation  with  respect  to  the 
matrix.  More  details  of  the  TEM  analyses  are  given  in 
the  Supporting  Information. 

Fig.  2d  shows  the  energy-dispersive  X-ray  spectroscopy 
(EDS)  results  obtained  from  the  matrix  and  the  precipitate 
measured  from  adjacent  regions  of  similar  thickness.  EDS 
microanalysis  of  the  sulfur-rich  precipitates  is  complicated 
by  the  fact  that  the  solitary  sulfur  Ky  emission  line 
(2.307  keV)  overlaps  with  that  of  the  Pb  Ma  line 
(2.342  keV).  Therefore,  the  intensity  of  the  precipitate 


spectrum  was  normalized  to  that  of  the  matrix  spectrum 
at  the  Pb  Ly  line  at  10.56  keV.  It  is  clear  that  the  precipi¬ 
tates  are  compositionally  distinct  from  the  matrix  with  no 
Te  spectrum  detected.  The  intensity  of  the  overlapped  Pb 
My  and  S  Ky  peaks  at  2.3  keV  is  marginally  greater  for 
the  precipitate  compared  with  the  matrix.  Although  precip¬ 
itates  generate  sulfur  X-rays,  the  very  high  mean  atomic 
number  of  the  precipitate  would  result  in  strong  absorption 
for  such  a  low-energy  X-ray.  The  mass  balance  and  elec¬ 
tron  diffraction  analyses  confirm  sulfur-rich  precipitates. 

Table  2  presents  the  room-temperature  Hall  carrier  con¬ 
centration  ( nH  =  rH/e.RH)  for  composite  samples.  The  hole 
carrier  concentration  for  all  samples  is  ~8  x  1019  cm-3. 
The  sodium  concentration  was  restricted  to  1  at.%  in  the 
current  study  due  to  the  low  reported  maximum  solubility 
of  sodium  in  PbTe  [32,33],  However,  the  room-temperature 
hole  carrier  concentration  for  the  samples  in  Table  2  is 
~8  x  1019  cm-3,  which  is  lower  than  the  optimum  carrier 


Table  2 

The  Hall  carrier  concentration  of  Na-doped  (PbTe)(o.9-*)(PbSe)o.i(PbS)x 
(.x  =  0.1,  0.15,  0.2,  0.25)  composite  samples  at  room  temperature. 


Sample 

Carrier  concentration  (cm  3) 

x  =  0.10 

8.3  x  1019 

x  =  0.15 

8.1  x  1019 

x  =  0.20 

7.3  x  1019 

x  =  0.25 

7.8  x  1019 
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Fig.  3.  Thermoelectric  transport  properties  of  (PbTe)(o.9_x)(PbSe)o.i(PbS)x  ( x  =  0,  0.05,  0.1,  0.15,  0.2,  0.25)  sintered  bulk  samples  doped  with  1  at.%  Na: 
(a)  electrical  resistivity  (mQcm),  cooling  curves;  (b)  Seebeck  coefficient  (pVK-1)  temperature  dependence;  (c)  measured  total  thermal  conductivity,  k 
(W  nr1  K_1),  and  calculated  lattice  thermal  conductivity  kl,  with  the  Lorenz  factor  approximated  using  a  single  parabolic  band  model  with  acoustic 
scattering;  (d)  electrical  resistivity  (mfl  cm)  of  Na-doped  (PbTe)co.9-JC)(PbSe)o.i(PbS)x  (x  =  0.1,  0.15,  0.2,  0.25)  sintered  bulk  sample  during  heating. 


concentrations  for  p-type  Pb-chalcogenides  [4,6,10,16,34], 
The  samples  are  highly  degenerate  semiconductors,  while 
enhanced  thermoelectric  efficiency  is  anticipated  for  the 
current  study  compounds  at  optimum  hole  carrier 
concentrations. 

The  electrical  resistivity,  Seebeck  coefficient,  measured 
total  thermal  conductivity,  k,  and  calculated  lattice  thermal 
conductivity,  kl,  for  sodium-doped  (PbTe)(o.9-X)(PbSe)0.i 
(PbS)Xj  0  =  0,  0.05,  0.1,  0.15,  0.2  and  0.25)  composite  sam¬ 
ples  as  a  function  of  temperature  in  the  range  of  300-850  K 
are  compared  with  single-phase  sodium-doped  (PbTe)(0.9_x) 
(PbSe)0.i(PbS)Xi  (x  =  0  and  0.05)  samples  [20]  in  Fig.  3.  All 
samples  show  behavior  typical  of  strongly  degenerate  semi¬ 
conductors,  with  the  Seebeck  coefficient  and  electrical  resis¬ 
tivity  increasing  with  temperature.  It  should  be  noted  that 
the  resistivity  curves  for  all  samples  containing  the  PbS- 
rich  secondary  phase  show  different  behavior  on  heating 
and  cooling  over  a  temperature  range  of  ~500-700  K. 
The  selected  heating  and  cooling  rates  are  100  °Ch-1. 
Fig.  3a  shows  the  smooth  decrease  in  resistivity  as  a  func¬ 
tion  of  temperature  for  Na-doped  (PbTe)(0.9_x)(PbSe)o.i 
(PbS)x>  (x  =  0.1,  0.15,  0.2  and  0.25)  samples  upon  cooling 
from  850  K.  The  resistivity  curves  of  single-phase  samples 
(x  =  0  and  0.05)  indicate  the  same  trend  during  cooling 


and  heating  [20],  The  resistivity  of  composite  samples, 
(PbTe)(0.9-X)(PbSe)o.i(PbS)x  (x  =  0.1,  0.15,  0.2  and  0.25), 
during  heating  in  Fig.  3d  show  relatively  distinct  behavior 
for  samples  with  a  distinct  secondary  phase  fraction.  The 
cycle  was  repeated  several  times,  and  various  heating  rates 
were  employed  to  confirm  the  presence  of  hysteresis  in  the 
resistivity  for  the  composite  samples.  Although  the  peak 
temperature  and  the  slope  of  the  curve  vary  slightly  with 
the  heating  rate,  the  cooling  curve  trends  were  independent 
of  the  rate.  This  behavior  might  be  attributed  to  the 
sodium  partitioning  [30,32],  and  the  changes  in  the  distri¬ 
bution  and  morphology  of  secondary  phases  [35],  This  is 
discussed  in  more  detail  elsewhere.[36] 

Fig.  3a  shows  that  the  electrical  resistivity  is  increased 
considerably  by  sulfide  addition  up  to  the  solubility  limit 
(x  <  0.1).  However,  sulfide  additions  above  this  level, 
which  increases  the  precipitate  fraction,  cause  only  minor 
further  increase  in  electrical  resistivity.  The  electrical  resis¬ 
tivity  of  single-phase  alloys  is  increased  by  sulfide  addition 
due  to  the  scattering  of  carriers  from  disordered  atoms  [37], 
and  possibly  due  to  the  increase  in  the  effective  mass  of  the 
light-band  by  alloying  with  sulfur. 

In  our  recent  report  [20],  we  indicated  that  the  addition 
of  PbS  to  sodium-doped  solid  solution  PbSe-PbTe  up  to 
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the  solubility  limit  provides  a  considerably  lower  thermal 
conductivity,  a  heavier  density  of  states  effective  mass,  a 
wider  band  gap,  and  a  larger  energy  offset  between  valence 
bands.  This  results  in  a  higher  temperature  for  band  con¬ 
vergence  and  a  larger  Seebeck  coefficient  in  sulfur-containing 
samples  and  provides  a  thermoelectric  efficiency  superior  to 
binary  PbQ  (Q  =  Te,  Se,  S)  and  ternary  systems  of 
PbTe-PbSe  and  PbTe  PbS.  These  results  have  been  con¬ 
firmed  by  a  more  recent  report  on  single-phase  p-type  qua¬ 
ternary  (PbTe)(  i  _2x/PbSe)dPbS)  Y  alloys  [21],  The  chemical 
compositions  of  the  matrix  and  the  second  phase  are  fixed 
at  any  given  temperature  for  nanocomposite  samples 
(x  >0.1)  and  the  weight  fraction  of  the  secondary  phase 
increases  with  x;  while,  as  the  temperature  increases,  the 
solubility  limit  of  PbS-rich  precipitates  increases  in  the 
matrix  of  composite  samples,  and  the  chemistry  of  material 
is  changed.  The  nanostructuring  through  transformation 
[17,38,39]  and/or  precipitation  [13,16]  generally  involves  a 
significant  fraction  of  secondary  phases.  The  physical  prop¬ 
erties  of  a  composite  material  are  related  to  the  volume 
fraction  and  the  physical  properties  of  the  individual  com¬ 
ponents  of  the  composite,  and  various  equations  have  been 
developed  to  model  these  associations  [40,41],  It  is  beyond 
the  scope  of  the  current  study  to  quantitatively  describe  the 
contribution  of  secondary  phases  on  electronic  transport 
properties  of  the  current  study  composites,  suffice  it  to  state 
that  the  chemical  composition  deviation  at  high  tempera¬ 
tures  along  with  variation  of  secondary  phase  fractions 
results  in  unpredicted  behavior  of  resistivity  and  Seebeck 
coefficient  of  composite  samples  with  a  large  fraction  of 
secondary  phases  (x  =  0.15,  0.20  and  0.25). 

The  thermal  conductivities  of  single-phase  quaternary 
Pb-chalcogenides  which  are  alloyed  marginally  with  PbS 
(up  to  10  at.%)  are  considerably  lower  than  parent  com¬ 
pounds  of  (PbTe)(i_x)(PbSe)x  [20,21],  Surprisingly,  the 
room-temperature  thermal  conductivity  of  nanostructured 
samples,  which  contains  secondary  sulfide  phase  (Fig.  3c), 
is  only  slightly  lower  than  the  thermal  conductivity  of  a  sin¬ 
gle-phase  alloy  of  (PbTe)o.85(PbSe)o.i(PbS)0.o5  which  con¬ 
tains  5  at.%  soluble  PbS.  This  difference  is  negligible  at 
higher  temperatures.  The  lower  thermal  conductivities  of 
the  phases  containing  secondary  sulfide,  compared  to 
PbS-free  single-phase  alloys  might  be  attributed  to  phonon 
scattering  on  solute  atoms  with  large  mass  contrast  [42] 
and/or  phonon  scattering  at  defects  and  interfaces  origi¬ 
nated  from  the  distributed  sulfide  secondary  phase  [16]. 
In  order  to  fully  realize  this  effect,  the  lattice  thermal  con¬ 
ductivity,  kl,  was  obtained  by  subtracting  the  electronic 
component,  Ke.  The  value  of  the  charge  carrier  thermal 
conductivity  Ke  can  be  determined  via  the  Wiedemann- 
Franz  relation,  Ke  =  LT/p,  where  p  is  the  resistivity,  and 
L  is  the  Lorenz  number  estimated  as  a  function  of  temper¬ 
ature,  assuming  a  parabolic  band  with  acoustic  phonon 
scattering  [43],  This  rough  estimation  has  been  shown  to 
be  reasonably  consistent  with  a  more  detailed  model  calcu¬ 
lation  taking  the  band  parabolicity  and  multiband  conduc¬ 
tion  effects  into  account  [4]  (see  Supporting  Information). 


Table  3 

The  total  and  lattice  thermal  conductivity  of  (PbTe)(o.9-x)(PbSe)o.i(PbSL 
( x  =  0,  0.05,  0.1,  0.15,  0.2,  0.25)  sintered  bulk  samples  doped  with  1  at.% 
Na  at  room  temperature  and  850  K. 


Sample 

Room-temperature 
lattice  thermal 
conductivity 

Room-temperature 
total  thermal 
conductivity 

Lattice  thermal 
conductivity  at 
850  K 

x  =  0 

1.4 

3.3 

0.7 

x  =  0.05 

1.26 

2.5 

0.6 

x  =  0.10 

1.2 

2.3 

0.6 

x  =  0.15 

1.2 

2.2 

0.6 

x  =  0.20 

1.4 

2.3 

0.6 

x  =  0.25 

1.5 

2.4 

0.6 

A  low  room-temperature  lattice  thermal  conductivity  of 
~1.2W  m-1  K-1  is  observed  for  the  sample  with  x  =  0.15, 
which  is  reduced  to  ~0.6  W  m-1  K-1  at  tempera¬ 
tures  >  650  K,  whereas  the  lattice  thermal  conductivity  is 
increased  for  samples  with  x  =  0.2  and  0.25.  Table  3  sum¬ 
marizes  the  total  and  lattice  thermal  conductivities  of  Na- 
doped  (PbTe)(0.9_X)(PbSe)o.i(PbS)x  (x  =  0,  0.05,  0.1,  0.15, 
0.2  and  0.25)  samples  at  room  temperature  and  850  K. 
The  phonon  scattering  on  solute  atoms  with  large  mass 
contrast  [42]  and  on  interfaces  of  nanoprecipitates  [16] 
reduces  the  lattice  thermal  conductivity  of  alloys  contain¬ 
ing  sulfur  up  to  x  =  0.15.  However,  further  addition  of  sul¬ 
fides  to  the  quaternary  system  results  in  a  larger  secondary 
phase  volume  fraction  in  the  composite.  The  lattice  thermal 
conductivity  of  PbS  is  higher  than  that  of  PbTe  [44],  and  a 
large  volume  fraction  of  PbS  precipitates  in  samples  with 
x  >  0.2  results  in  a  subsequent  increase  in  their  lattice  ther¬ 
mal  conductivity. 

Fig.  4a  compares  measured  zT  up  to  850  K  for  compos¬ 
ite  samples  of  (PbT e)(()  9_ x)(PbSe)0.  i (PbS)x  (x  =  0.1,  0.15, 
0.2  and  0.25)  to  single-phase  alloys  with  x  =  0  and  0.05 
[20].  The  maximum  zT  value  of  1.6  was  achieved  at 
750  K  for  the  single-phase  sample  with  x  =  0.05,  while 
the  thermoelectric  efficiency  is  reduced  for  composite  com¬ 
pounds.  Although  a  high  thermoelectric  efficiency  of  ~1.4 
has  been  achieved  for  the  composite  sample  with 
x  =  0.25,  the  figure  of  merit  of  composite  compounds 
and  alloys  with  a  secondary  phase  are  lower  for  than  sin¬ 
gle-phase  alloys.  This  efficiency  is  obtained  in  the  sample 
where  35%  of  Te  atom  sites  are  occupied  by  the  abundant 
elements  selenium  and  sulfur. 

The  thermoelectric  figure  of  merit  of  the  p-type  compos¬ 
ite  (PbTe)o.65(PbSe)o.i(PbS)o.25  in  the  current  study  is 
compared  with  the  maximum  zT  values  reported  for  p-type 
Te-free  PbSe  [10],  Sr-added  PbS  [12],  CdS-added  PbSe  [45], 
single-phase  PbTe  [34],  Te-reduced  PbTe-12  at.%  PbS  [16] 
and  PbTe-16at.%  PbS  [16]  composites  with  identical 
dopant  concentration  in  Fig.  4b.  A  figure  of  merit  of  >1 
has  been  obtained  at  temperatures  >  600  K,  compared  with 
750  K  for  p-type  PbSe  [10]  and  800  K  for  p- type  Sr-added 
PbS  [12],  The  high  zT  obtained  at  lower  temperatures  can 
enhance  the  power  generation  efficiency  of  thermoelectric 
devices.  The  nanocomposite  sample  containing  only 
32.5  at.%  Te,  (PbTe)0.65(PbSe)0.i(PbS)0.25,  shows  similar 
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Fig.  4.  (a)  Temperature  dependence  of  the  dimensionless  figure  of  merit  zT  for  1  mol.%  Na-doped  (PbTe)(0.9_x)(PbSe)o.i(PbS)JC  ( x  =  0,  0.05,  0.1,  0.15,  0.2, 
0.25)  sintered  bulk  samples;  (b)  comparison  of  the  thermoelectric  efficiency  of  the  nanocomposite  Pbo.98Nao.02S0  25Seo.1Teo.65  quaternary  system  with 
maximum  reported  values  for  binary  p-type  PbSe  [10],p-type  Sr-added  PbS  [12],p-type  CdS-added  PbSe  [45],  and p-type  PbTe  [34],  PbTe-12  at.%  PbS  [16] 
and  PbTe-16  at.%  PbS  [16],  indicating  larger  efficiency  over  a  wide  temperature  range. 


efficiency  to  PbTe  [34]  and  (PbTe)0.88(PbS)0.i2  [16]  with 
much  higher  tellurium  content,  while  it  has  a  considerably 
higher  thermoelectric  efficiency  than  PbSe  [10,45],  Sr-added 
PbS  [12]  and  (PbTe)0.84(PbS)0.i6  [16]  with  lower  tellurium 
content. 

4.  Conclusion 

In  summary,  we  demonstrate  high  thermoelectric  effi¬ 
ciency  in  p- type  composite  quaternary  Pb-chalcogenides 
where  tellurium  is  replaced  by  the  abundant  elements  sele¬ 
nium  and  sulfur.  The  partial  dissolution  of  sulfide  precipi¬ 
tates  in  the  matrix  at  elevated  temperatures  results  in  a 
higher  Seebeck  coefficient  than  in  single-phase  alloys  in 
combination  with  reasonably  low  electrical  resistivity. 
Our  results  show  that  the  very  low  achieved  thermal  con¬ 
ductivity  is  dominated  by  phonon  scattering  at  high  con¬ 
trast  atomic  mass  solute  atoms  rather  than  at  distributed 
secondary  phases  within  the  matrix.  We  show  that  the  ther¬ 
mal  conductivity  of  composites  with  a  considerable  frac¬ 
tion  of  sulfide  precipitates  is  higher  than  for  single-phase 
sulfur-alloyed  compounds,  but  still  lower  than  for  the  sul¬ 
fur-free  compound.  Our  results  are  a  step  forward  in  real¬ 
izing  the  thermoelectric  performance  of  multiphase  alloys 
in  quaternary  Pb-chalcogenides  with  lower  tellurium 
concentration. 
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